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Quantum simulation promises to be one of the primary application of quantum computers, should 
one be constructed. This article briefly summarizes the history quantum simulation in light of the 
recent result of Wang and coworkers demonstrating calculation of the ground and excited states for 
a HeH+ molecule, and concludes with a discussion of why this and other recent progress in the held 
suggests that quantum simulation of quantum chemistry has a bright future. 


In 1982, the physicist Richard Feynman noted that it is 
very hard for a classical computer to simulate some quan¬ 
tum mechanical systems - but that if “the computer it¬ 
self be built of quantum mechanical elements which obey 
quantum mechanical laws” [l[ , then the simulation would 
become far easier. With this suggestion, Feynman inau¬ 
gurated the field of quantum computing, and introduced 
the idea that quantum computers change the boundaries 
of computational complexity. 

Simulation of quantum physical systems is thus the 
first “killer app” of quantum computing. Since then, 
several other fast quantum algorithms have emerged, no¬ 
tably Peter Shor’s algorithm for integer factorization and 
discrete logarithms 13 and Lov Grover’s algorithm for un¬ 
structured search (3|. But quantum simulation remains 
one of the primary applications of a quantum computer, 
should one be constructed. It is arguably the most in¬ 
triguing and potentially valuable of all the known quan¬ 
tum algorithms, with potential impact in chemistry 4, 
materials science Q , and elementary particle physics ^ ■ 

The canonical quantum simulation algorithm uses the 
same subroutine as Shor’s integer factorization algorithm 
- the quantum Fourier transform - but instead of us¬ 
ing it to reveal the period of a modular function, sim¬ 
ulation uses it to estimate the energies of a Hamilto¬ 
nian. It does this using Kitaev’s algorithm for quantum 
phase estimation [^, which efficiently reveals the phase 
(and thus the eigenvalue), of an eigenvector of a uni¬ 
tary operation. Abrams and Lloyd [ 3 , Q showed that a 
broad class of Hamiltonians could be efficiently simulated 
(and thus probed using phase estimation) using a Suzuki- 
Trotter expansion. They applied their results to the Hub¬ 
bard model and suggested that these techniques would 
lead to an exponential speed-up over classical comput¬ 
ing resources. Somma and coworkers Q extended quan¬ 
tum simulation to many more systems by showing that 
the Jordan-Wigner transformation 0 could be used to 
map the creation/annihilation operators that often de¬ 
fine quantum mechanical Hamiltonians into N-qubit op¬ 
erators. This paved the way for quantum simulation of 
quantum chemistry. 

Quantum chemistry addresses the problem of describ¬ 
ing the electronic structure of molecules and materi- 
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als, which provides information about how they dis¬ 
sociate, react, absorb light, and interact with other 
molecules and materials. Typically, quantum chemical 
methods describe molecular electronic states as prod¬ 
ucts of single-electron states generated from a variety of 
effective Hamiltonians, including Hartree-Fock and 
density functional theory [T^ . [l3l | . A huge variety of in¬ 
teresting chemical systems are well-described either by 
these product states or by simple perturbations of them. 
However, these methods break down in systems involving 
strongly correlated electrons. For these molecules, quau- 
tum chemists fall back ou brute force methods, the most 
accurate (and arduous) of which is full configuration in¬ 
teraction (Cl) [l^ . 

Full Cl considers the many-body Hamiltonian formed 
by all the possible ways of distributing N electrons into 
M one-electron states. Because the number of such con¬ 
figurations grows combinatorially with N and M, full Cl 
computatious require [classical] computational resources 
that scale exponentially with the molecule’s size. This 
kind of scaling rapidly surpasses the capacity of auy exist- 
iug computer, aud so quantum chemists have developed 
a variety of “truncated” approximations to full CL Typ¬ 
ically, these involve configurations that can be generated 
using just a few excitations from the ground state elec¬ 
tronic configuration. Unfortunately, although truncated 
Cl wave functions can be evaluated with only polynomial 
resources, they are known to lack properties, such as size 
consisteucy [l5| , that are required for quantitative chemi¬ 
cal predictions. Related perturbational approaches, such 
as coupled-cluster approaches [l^, iuclude size consis¬ 
tency, but often lack strict variational bounds on the en¬ 
ergy. So, although both truncated Cl and coupled cluster 
methods often give excellent energies in practice, they are 
incomplete solutions to chemical simulation because they 
lack certain properties of the full Cl method. For the 
chemical problems where these methods fail, fully quau- 
tum simulation a la Feynman might be the only viable 
approach. 

In 2005, Aspuru-Guzik and coworkers applied iter¬ 
ative phase estimation and other techniques from quan¬ 
tum information theory to the full Cl problem. Although 
their approach scales polynomially for the full Cl prob¬ 
lem, the exact scaling of their algorihtm was unclear in 
2005. However, their approach suggested that even mod¬ 
est quantum resources could poteutially euable full Cl 
simulations that outperform the largest, fastest super- 
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computers. 

The first experimental demonstration of quantum full 
Cl algorithm came 5 years later, when a photonic quan¬ 
tum information processor [l^ calculated the ground 
and electronically excited states for a minimal basis set 
description of H 2 using 20-bit iterative phase estima¬ 
tion. Although this proof-of-principle computation did 
not provide any added insight into the nature of the 
H 2 bond, it demonstrated conclusively that the quantum 
phase estimation algorithm could produce accurate bond¬ 
ing and excitation energies - even in the face of the noise 
and decoherence inherent to imperfect physical qubits. 

In this issue of ACS Nano., Wang and coworkers [l^ 
demonstrate a quantum calculation of the ground and ex¬ 
cited states for a HeH+ molecule, which is isoelectronic 
with H 2 , but has differently-charged nuclei. Wang et al. 
report the highest precision achieved to date in the quan¬ 
tum simulation of molecular energies; they surpass chem¬ 
ical precision by 10 orders of magnitude. Their quantum 
simulation was performed using two qubits from a dia¬ 
mond NV center, and represents the first implementation 
of the full Cl algorithm on a solid-state qubit. As with 
the earlier demonstration [I^ , this result is most notable 
for verifying that these algorithms can be successfully 
implemented on actual qubits, rather than just in theo¬ 
retical idealizations. 

The principle of quantum chemistry simulation has 
been proved. The next urgent question is how it will 
scale - how many qubits and how much time will be re¬ 
quired to apply quantum simulation algorithms to larger, 
more general molecules. The time requires scales with iV, 
the number of basis functions used to describe the atomic 
orbitals that comprise the molecule (adequate basis sets 
require 5-20 basis functions per atom), and/or Z, the 
maximum nuclear charge in the molecule. The initial 
analysis last year suggested that quantum simulation 
algorithms might require 0{N^) clock cycles. Although 
this is a huge advance, in-principle, over the 0{e^) scal¬ 
ing of exact classical simulation, it remains prohibitively 
intractable in practice. However, more careful analyses of 
the Suzuki-Trotter expansion have produced steady and 
rapid improvements in this scaling, first to 0{N'^) [2l|, 
then to [ 2 ^, and, most recently, to 0{N^Z^-^) 

M- 

The best scaling known today suggests that, given a 
handful of good qubits, quantum simulation could pro¬ 
duce practical, useful results. A closer look, however, 
suggests a less optimistic picture in the near future. 
Molecules with N = 1000 are routinely analyzed on lap¬ 
top computers using density functional theory (DFT), 
the workhorse method of chemistry and materials sci¬ 
ence. The most optimistic scaling given above suggests 


that a quantum computer would need at least 10® op¬ 
erations to match what DFT can do on a laptop. Be¬ 
cause qubits have far higher error rates than classical 
computers, such a large computation would absolutely 
demand quantum error correction Error correction 
imposes massive overhead in time, number of qubits, and 
complexity, because small quantum rotations have to be 
compiled into H, S, and T gates [ 2 ^, some of which must 
be exhaustively distilled from noisy resources [ 2 ^ . These 
considerations suggest that practical quantum simulation 
of meaningful molecules is substantially more challenging 
than it appears at first - and perhaps infeasible in the 
near term. 

Perhaps surprisingly, we remain optimistic. Although 
the demonstrations by Wang et al [T^ and others are a 
long way from practical utility, they demonstrate that 
the principle is sound. It is widely believed that fur¬ 
ther algorithmic improvements are possible - after all, 
just in the past two years we have seen the time scal¬ 
ing drop from 0{N^) to 0{N^Z^-^). Moderate improve¬ 
ments in scaling could dramatically change feasibility, 
and alg orithms that avoid the Suzuki-Trotter expansion 
[27H30| appear promising. Finally, we observe that some 
of the most promising candidates for near-term “quan¬ 
tum supremacy” (practical speedups over existing clas¬ 
sical computers) are robust quantum simulation algo¬ 
rithms that may not require error correction, either by 
using device noise to mimic real-world noise in the simu¬ 
lated system, or by leveraging shallow quantum circuits 
that run so quickly that errors do not accumulate. Ana¬ 
log quantum simulation is one such candidate, and has 
been used to mimic Hubbard 311 and spin [32j| Hamil¬ 
tonians. Analog methods are not currently known for 
chemical or molecular systems, but these systems can 
be addressed within the variational eigensolver approach 
[ 3 ^, which leverages a small quantum computer to eval¬ 
uate the energy terms of a parametric wave function that 
is varied by an associated classical computer-driven opti¬ 
mization. We believe that between the near-term promise 
of these non-traditional algorithms, and steady progress 
toward the long-term goal of digital quantum simulation, 
quantum simulation of quantum chemistry has a bright 
future. 
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